Alteration of gas phase ion polarizabilities upon hydration in high dielectric liquids 
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We investigate the modification of gas phase ion polarizabilities upon solvation in polar solvents 
and ionic liquids. To this aim, we develop a classical electrostatic theory of charged liquids composed 
of solvent molecules modeled as finite size dipoles, and embedding polarizable ions that consist of 
Drude oscillators. In qualitative agreement with ab-initio calculations of polar solvents and ionic 
liquids, the hydration energy of a polarizable ion in both type of dielectric liquid is shown to favor 
the expansion of its electronic cloud. Namely, the ion carrying no dipole moment in the gas phase 
acquires a dipole moment in the liquid environment, but its electron cloud also reaches an enhanced 
rigidity. We find that the overall effect is an increase of the gas phase polarizability upon hydration. 
In the specific case of ionic liquids, it is shown that this hydration process is driven by a collective 
solvation mechanism where the dipole moment of a polarizable ion induced by its interaction with 
surrounding ions self-consistently adds to the polarization of the liquid, thereby amplifying the 
dielectric permittivity of the medium in a substantial way. We propose this self-consistent hydration 
as the underlying mechanism behind the high dielectric permittivities of ionic liquids composed of 
small charges with negligible gas phase dipole moment. Hydration being a correlation effect, the 
emerging picture indicates that electrostatic correlations cannot be neglected in polarizable liquids. 

PACS numbers: 05.20.Jj,61.20.Qg,77.22.-d 



I. INTRODUCTION 

The atomic electron cloud distortion induced by an 
external field is strongly influenced by the dielectric en- 
vironment embedding the atom. This distortion ability 
referred as the induced polarizability is one of the key 
ion specific effects in the simulation of salt solutions in 
inhomogeneous media such as the water-air interface or 
protein-water surfaces p_J . The precise knowledge of the 
change in the polarizability of an isolated ion upon hydra- 
tion in water is particularly important for the develop- 
ment of polarizable force fields used in these simulations. 
Moreover, ionic polarizability is also believed to have a 
substantial effect on the polarity of ionic liquids. Indeed, 
numerical studies based on ab-initio calculations show 
that the large dielectric permittivity of ionic liquids such 
as [C2mim] [NTf 2] and [C2mim][BF4] composed of ions 
with small individual dipole moments cannot be solely 
explained by their rotational polarizability [2]. This sug- 
gests that an additional polarization mechanism result- 
ing from the interaction of the polarizable ion with the 
surrounding ions in the liquid must be present. 

The alteration of ionic gas phase polarization upon sol- 
vation has been so far considered within numerical ap- 
proaches based on quantum calculations with polarizable 
continuum model (PCM) or explicit solvent. These two 
approaches interestingly yield diverging pictures on the 
hydration of polarizable ions. Namely, the calculations 
with explicit solvent indicate that the ionic polarizability 
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is decreased with respect to the gas phase [3], whereas 
PCM approaches yield a higher polarizability in the liq- 
uid state [U [5] (see also Ref. [B] for a review on the com- 
putational state of the art). The latter case is also in line 
with the ab-initio calculations of pure water clusters [7] 
and ionic liquids (8J, where the transfer of both type of 
molecules from gas to the liquid environment was shown 
to increase their dipole moment. 

In order to understand the physics behind the hydra- 
tion of polarizable molecules, analytical theories offer- 
ing a deeper understanding are needed. The theoreti- 
cal formulation of the problem requires in turn an ex- 
plicit and realistic consideration of the discrete charge 
structure of solvent molecules and ions. Unfortunately, 
this level of refinement has been until recently beyond 
the state of the art of electrostatic theories, which are 
mostly based on dielectric continuum solvents embedding 
point charges. The first statistical theory of inhomoge- 
neous electrolytes with explicit solvent was introduced in 
Ref. [5] in the form of a mean-field (MF) dipolar Poisson- 
Boltzmann (DPB) equation. This approach that models 
the solvent molecules as point dipoles was later gener- 
alized by including the steric interactions between the 
particles for inhomogeneous charged liquids [TO], and a 
one-loop extension was presented as well in Ref. [IT] to 
explain the salt induced dielectric decrement effect in 
bulk electrolytes. We have recently incorporated into 
the DPB approach surface polarization effects, which al- 
lowed us to significantly improve the agreement of the 
dielectric continuum electrostatic with experimental ca- 
pacitance data of carbon based materials [12 J. Sophisti- 
cated electrostatic formulations accounting for the dipo- 
lar and higher order multipolar moments of ions in the 
point dipole limit have been also proposed in Refs. [TBT - 



16] . In a similar context, we can also mention the works 
of Refs. p"7HT9] where the extended charge structure of 
rigid linear molecules was ingeniously considered. 

We have recently developed a non-local electrostatic 
theory of polar liquids with explicit solvent and polar- 
izable ions beyond the point dipole approximation |20j . 
The electrolyte model that treats solvent molecules as 
finite size dipoles and polarizable ions as Drude oscil- 
lators was investigated at the MF level. It was shown 
that the consideration of the extended charge structure 
of solvent molecules enables us to capture the non-local 
dielectric response of water at charged interfaces observed 
in molecular dynamics simulations and atomic force ex- 
periments. In this article, we reconsider the model of 
Ref. [20] beyond the MF level of approximation in or- 
der to characterize the hydration induced modification of 
ionic polarizabilities in high dielectric bulk liquids. We 
review in Section |TT] the derivation of the field theoretic 
charged liquid model, and derive the closure equations 
accounting for the correlations between the ions and the 
solvent molecules. These equations are first solved in 
Section |III A| in order to investigate the hydration of 
a single polarizable ion in a polar solvent such as wa- 
ter. Then, within the same theoretical framework, we 
consider in Section |IHB| an ionic liquid free of solvent 
molecules in order to investigate a collective polarization 
effect in the liquid. It is shown that in both systems, 
our simple theory can capture the solvation induced elec- 
tronic cloud expansion effect observed in ab-initio calcu- 
lations [H 13 |H] , and provides a physical explanation in 
terms of the electrostatic energy released by the ion upon 
hydration. The limitations of the liquid model and the 
computation scheme, and necessary extensions are dis- 
cussed in detail in the Conclusion. 



Drude oscillator model 



II. MODEL 



We briefly review in this section the derivation of 
the field theoretic partition function for the polar liquid 
model previously introduced in Ref. [20 . Then, starting 
from the Dyson equation, we derive an integral equa- 
tion for the dielectric permittivity function embodying 
the interactions between the polarizable ions and solvent 
molecules of the bulk liquid. 

The geometry of solvent molecules is depicted in 
Fig.nTa). The polar liquid is composed of overall neutral 
solvent molecules modeled as linear dipoles of length a, 
and two point charges of valency ±Q = ±1 at the ex- 
tremities. Furthermore, the solvent contains polarizable 
molecules of p species, each of them being an oscillating 
rod of length b (see Fig. [IJb)). The point charges e^ and 
Cj at the extremities satisfy the inequality e^c,; < 0, where 
the index i — l...p runs over the ionic species. Moreover, 
the ionic polarizability is taken into account within the 
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where the square of the variance of electronic cloud oscil- 
lations bp t is proportional to the induced polarizability of 
ions a in the gas phase [20]. Because the former offers a 
more intuitive realization of the electronic cloud fluctua- 
tions induced by thermal excitations, we will discuss the 
results in terms of the length scale b P i. Furthermore, in 
the present work, we will consider exclusively the case of 
equal ionic polarizabilities for all species, but the analyti- 
cal results will be given for the general case. We also note 
that the electroneutrality condition implies the equality 
J2i PibQi = 0: with pa the bulk density, and qi — e, + c t 
the total charge of the polarizable molecules with species 
i. 

The canonical partition function for the system com- 
posed of solvent molecules and ions coupled with electro- 
static interactions read 
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where N s is the total number of solvent molecules, N^ the 
number of ions for the species i, and Xrd and \t{ denote 
respectively the thermal wavelengths of solvent molecules 
and ions. We also introduced in Eq. pi) the compact no- 
tation v = ({xfc}, {afe}, {yij}, {bj}) for the vectors char- 
acterizing the configuration of particles, with x^ and y^ 
denoting respectively the coordinate of the charges +Q 
and &i of the solvent molecules and polarizable ions in de- 
picted in Fig. [I] Furthermore, flk = (0k,fk) is the solid 
angle characterizing the orientation of the kth solvent 
molecule, 9 being the angle between the oriented dipole 
and the z-axis (see Fig. flTa)). We finally note that in 
Eq. (pi, we subtracted from the total Hamiltonian the 
self energies of ions and polar molecules in the air, e^ = 
(e| + cf ) u c (0)/2 + e lCl v c {b) and e s = Q 2 [v c (0) - v c (a)}. 
This point will be discussed below in further detail. 

The Hamiltonian of the bulk liquid is composed of pair- 
wise electrostatic interactions, 



H e i(v) = - / [p lc + Psc] r v c (r - r') [p ic + p sc ] r , , (3) 

where the total ionic and solvent density operators for 
the charge compositions depicted in Fig. [T] are defined as 

P N t 

ft c (r) = X!X! [ e ' ;(5 ( r ~ y ^') + Cl(5 ( r ~ Vi i - b J")]( 4 ) 
i=\ j=i 

p sc (r) = ^Q [6(r - Xfc) - <J(r - x fc - a*.)] . (5) 
fe=i 





FIG. 1: (Color online) Charge composition of solvent 
molecules of size a (a) and polarizable ions with a fluctuating 
length b (b). In the present work, we consider exclusively the 
case of ionic valencies a and a of opposite sign (ejCj < 0), 
and solvent molecules with monovalent point charges Q — 1. 



Moreover, in Eq. pi, v c (r — r') = £s/|r — r'| stands for 
the Coulomb potential in the air medium, with £g = 
e 2 / {AiteqUbT] ~ 55 nm the Bjerrum length and £o the 
dielectric permittivity in the air, e the electron charge, 
and T — 300 K the ambient temperature. We note that 
in the rest of the article, dielectric permittivities will be 
expressed in units of the air permittivity £rj, and energies 
in units of the thermal energy k B T. 

In order to transform the partition function (pi) 
into a more tractable form, we pass from the parti- 
cle density to the fluctuation potential representation 
by performing a standard Hubbard-Stratonovich trans- 
formation. In this representation, the grand canoni- 
cal partition function of the system defined as Zq = 

^N B >oHi=i^N i >o efliNieflwNs Z c takes tne f° rm of a 
functional integral - over the fluctuating electrostatic po- 
tential </>(r), Z G = fVcj) e~ H[4> \ with the Hamiltonian 
functional 20 
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The first term on the r.h.s. of Eq. p]) is the electro- 
static energy of the freely propagating field in the air. 
The second term corresponds to the density of solvent 
molecules, and their fugacity is denoted by A s . Finally, 
the third term on the r.h.s. of Eq. (pi) is the density of 
polarizable ions with fugacity A^. 

The Hamiltonian pi was already derived in Ref. [20] 
for the more general case of multipolar solvents embed- 
ding polarizable ions, and the saddle-point solution of 
the partition function corresponding to the MF approxi- 



mation was investigated for polar liquids in contact with 
charged planes. In order to account for correlation effects 
in the bulk liquid beyond the MF level, we need to derive 
the electrostatic correlation function. Our starting point 
is the following form of the Dyson equation, 
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where J(r) is a generalized current introduced for the 
derivation of the two point correlation function. A 
proof of the equality ^ can be found in Ref. [22] . We 
also remind that the derivation of the electrostatic self- 
consistent equations of the primitive ion model [23] with 
the use of this equality was presented in Ref. [21] • By tak- 
ing now the functional derivative of Eq. pi) with respect 
to J(r') and setting J(r') = 0, one obtains the following 
equation for the two point correlation function, 
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where the bracket (•) denotes the field average with the 
Hamiltonian Functional in Eq. pi) . In Eq. p| , the depen- 
dence of the fluctuating solvent and ion densities (i.e. the 
functions inside the brackets on the l.h.s.) on the values 
of the potential at different points around r is a signature 
of non-local electrostatic interactions resulting from the 
extended charge structure of the solvent molecules and 
ions [2D]. 

We emphasize that the formal equation pi is an exact 
relation. However, because the Hamiltonian of Eq. pi) 
is non-linear in the potential </>(r), an exact analytical 
evaluation of the averages over the fluctuating poten- 
tial is impossible. To progress further, we approximate 
this non-linear Hamiltonian with a quadratic Hamilto- 
nian functional, 
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where the electrostatic potential is chosen as the solution 
of the equation ^, that is, w (r,r') = (</>(r)(/>(r')) . At 
this stage, we note that the spherical symmetry in the 
bulk liquid implies t>o(r,r') = t>o(r — r'), and this allows 
us to expand the potential in Fourier space as i>o(r — r') 



/ 
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^ e 
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v a (q). Evaluating the averages in Eq. (8) 



with the quadratic functional (19]) and injecting into the 
result the Fourier expansion of the correlation function, 



the explicit form of the potential finally follows in the 
form [25] 
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with the Fourier transformed dielectric permittivity func- 
tion 
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We introduced in Eq. (|11|) the solvent and ionic screen- 
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ing parameters in the air, 

8ir\eiCi\£BPib- Furthermore, we defined in Eq. (Ill the 



statistical average over the fluctuations of the electronic 
cloud, 
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with the potential of mean force (PMF) 
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where the Fourier transform of the Coulomb potential in 
the vacuum given by v c (q) — q 2 /(AttIb)- We also note 
that deriving Eq. (10), we used the thermodynamic rela- 
tions between the particle fugacities and concentrations, 
p sb = A s d\nZ/dA s = A s e"^ and p ib = A.dlnZ/dAi = 
A^dhe-^^-^^/^b^f/ 2 , with the liquid state 
self-energies of solvent molecules and ions respectively 
defined as 



ips = - ( 
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One can notice that the energies in Eq. ( 14 1 and ( 15 ) cor- 



respond to the hydration energies of the solvent molecules 
and polarizable ions, i.e. the electrostatic cost to drive 
the molecules from the gas to the liquid environment. 
Moreover, one sees that Eqs. (131 and ( 15 ) are related as 
4>i{b) = —q 2 [v c (0) — v (0)] /2 + 4>ib(b), which indicates 
that the PMF ipip(b) brings the net contribution from 
the polarizability to the ionic hydration energy. Finally, 
unlike previous point dipole models where the electro- 
static energies have to be regularized with an ultraviolet 
cut-off in Fourier space [Til US] > our consideration of the 
finite solvent molecular size and electronic cloud exten- 
sion resulted in a cut-off free theory with well defined self 



energies in Eqs. (13)-(15) 



At this stage, we note that our motivation for sub- 
tracting from the Hamiltonian the gas phase self-energy 
of polarizable ions in Eq. pi) was twofold. First of all, this 
step allowed us to avoid the dipolar catastrophy problem. 
Indeed, the classical Drude oscillator model of Eq. (fll) 
does not prevent the electron from falling into the nu- 
cleus, and this results in divergent ionic self-energies for 
b — > 0. This problem could be avoided in an alternative 
way by modifying the Drude model with a cut-off at small 
b, but we found that this technical complication shadows 
the transparency of the analytical results. Furthermore, 
the Drude potential is clearly an approximative fashion 
to consider the quantum mechanical interatomic interac- 
tions that already include the electrostatic coupling be- 
tween the electron and the nucleus. We also note that 
the subtracted self-energy of solvent molecules does not 
affect the statistical average in Eq. (12). 



The relations ( 10 )-( 13 1 form a set of closure equations 
that should be solved self-consistently. These two rela- 
tions can be also interpreted as a single integral equation 
for the dielectric permittivity function e(q) in Fourier 
space. Then, one notes that computing the average in 



Eq. ( 11 1 by neglecting the PMF ( 13 1 in Eq. ( 12 ), one ob- 



tains the MF permittivity function derived in Ref. [2"U] . 
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Hence, electrostatic correlation effects are incorporated 
in the hydration PMF ipib(b)- In the rest of the article, 



the solution of the closure equations ( 10 )-( 13 1 will be con 



sidered in order to investigate the solvation of polarizable 
ions in high dielectric liquids. 



III. RESULTS 



In this section, we solve the closure equations ( 10 ) 



(131 in order to shed light on the electrostatic mecha- 



nism behind the hydration effects observed in ab-initio 
calculations for polarizable ions in high dielectric liquids 
such as polar solvents [H [5] and ionic liquids [21 18] . We 
first investigate in Section [III A| the hydration of a single 
polarizable ion in a polar liquid such as water, and we 



characterize in Section III B a similar cooperative solva- 
tion mechanism in ionic liquids exclusively composed of 
polarizable ions. 



A. Hydration of a single polarizable ion in water 

This section is devoted to the hydration of a single po- 
larizable ion in a strongly polar liquid such as water. In 
the dilute ion regime, the PMF of Eq. (13) has to be 



evaluated at the leading order in the ion concentration 
by neglecting the ionic contributions corresponding re- 
spectively to the second and third terms on the r.h.s. of 




FIG. 2: (Color online) Drude oscillator potential Eq. (fTl) (blue 
curve), hydration energy of Eq. ( 17 1 (red curve), and total dis- 



tortion potential of an hydrated polarizable ion (black curve) 
Model parameters are a — 1 A, p s t = 10 -4 M 



b P i =1A, and 



value tpib = — |ejCj| (re a a) £s/(2a). Thus, for dilute sol- 
vents, the hydration modifies the electronic cloud rigidity 
mainly at separation distances below the solvent molec- 
ular size. 

To extend the investigation of the hydration induced 
modification of the electronic cloud radius and rigidity 
beyond the dilute solvent regime, we can map Eqs. (10) 
and ( 13 ) onto an effective polarizable ion model. By ad- 



sorbing the effect of the hydration potential ipib(b) into 
an effective Drude oscillator model 



h(b) - (b ~ b "" )S 



(18) 



with the average dipole moment (or electronic cloud ra- 
dius) b mi and induced ion polarizability b V i in the liq- 



uid environment, and evaluating the average in Eq. ( 11 ) 



with the distortion potential ( 18 ) without the hydration 



PMF ( 13 1, we are left with the effective permittivity func- 
tion 



Eqs. (10) and (11). In order to illustrate the hydration 



mechanism in an intuitive way, we first consider a polar- 



izable ion in a dilute solvent. By expanding Eq. (13) at 
the order O ((ft s a) 
ular size a 



2 ^ , which is valid for the solvent molec- 
1 A in the solvent density regime p s b < 0.1 



M, one obtains for the PMF associated with the polariz- 
ability the close form expression, 
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The hydration potential ipi P (b) of Eq. (17) and the to- 
tal distortion energy hi(b)+tpi P (b) are compared in Fig. [2] 
with the distortion potential of an isolated ion hi (b) . One 
sees that the negative hydration potential ipi P {b) results 
in a net reduction of the bare distortion energy hi (b) . In 
other words, the hydration of a polarizable ion favors the 
expansion of its electronic cloud. This peculiarity results 
from the fact that the Born energy of a point charge is 
proportional to the square of its valency, and the point 
charges on the polarizable ion are of opposite sign and 
satisfy the inequality ef + c 2 > q 2 . As a result, the sol- 
vation energy of two separate charges with valencies e^ 
and Ci is lower than the Born energy of a single ion of 
valency qi in Eq. (15), that is ipi(b —> oo) < ipi(b = 0). It 



follows from this remark that for a rodlike molecule with 
the charges e^ and Ci of the same sign, hydration would in 
turn lead to a compression of the electronic cloud. Fur- 
thermore, the black curve in Fig. [2] shows that the total 
distortion potential hi(b) + ipip(b) exhibits a minimum. 
This means that the polarizable molecule without aver- 
age dipole in the gas phase acquires a net dipole moment 
upon hydration. One finally notes that in Eq. (17), the 



hydration potential converges for b > a to a constant 
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The comparison of the function ( 19 1 with Eq. ( 16 1 indi- 



cates that at the MF level, the ion has no dipole moment 
{bmi = 0), and its polarizability is equal to the gas phase 
value (b v i = b p i). By expanding now Eqs. ( 11 ) and ( |19| in 
the infrared (IR) regime up to the order G(q 4 ) and iden- 
tifying the quadratic and quartic terms in the wavevector 
q, one obtains the coupled equations 6b 2 vi + b mi = <V) 
and 60b 4 m + 20b 2 mi b 2 m + b A ml = (b 4 ). The solution of these 
equations respectively yields for the average dipole mo- 
ment and induced polarizability of the hydrated ion 
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where we introduced the total ionic polarizability 
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We evaluated the dipole moment and polarizabilities in 
Eqs. (20)-(22) with the numerical solution of Eqs. (10)- 
( 13 1 . Figure|3Ja) displays the variation of the ionic dipole 



moment b m i with solvent density for the gas phase po- 
larizability b p i — 0.2 A and various molecular valencies 
(solid curves). First of all, it is seen that an increase of 
the solvent concentration is accompanied with a mono- 
tonic rise of the dipole moment from zero to b m i ~ 4 A, 
until the latter saturates in the density regime p s b > 10 M 
where the ion becomes fully hydrated. Then, one notices 
in Fig. |3tb) that the expansion of the average electronic 
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FIG. 3: (Color online) (a) Enhancement of the ionic dipole moment b m i introduced in Eq. (20 1 and (b) the total polarizability 

btot,i 

density. The ion in t 



defined in Eq. (22 1, and (c) the reduction of the effective intrinsic polarizability b v i of Eq. ( 21 ) with increasing solvent 

0.20 A, and different ionic valencies from leicA — 1 to 



le polar solvent has gas phase polarizability b p 



3 are considered. The results obtained from the numerical solution of the self-consistent equations (101-1131 at a fixed dipole 
moment po — 1 A are displayed by solid curves for the solvent molecular size a = 1 A and ion valencies |ejC»| = 1 to 3, and by 
dash-dotted black curves for a = 3 A and divalent molecules with |ejCj| = 3. Dotted curves in (a) and (b) denote for divalent 
molecules the point dipole results of Eq. ( |35[l obtained in the limit a — > of Eqs. ( 10 1- ( 13 1 at fixed dipole moment, and circles 
mark the asymptotic equations (36 1 and (|37[) derived in the same point dipole limit for large solvent concentrations. Dashed 



horizontal curves correspond to the complete ionic hydration state of Eqs. (28 1- (30 1, (d) Dielectric permittivity profile around 
a point ion at r = for the solvent density p s t = 55 M. 



cloud radius upon hydration results in turn in an ampli- 
fication of the total polarizability b to t,i by several factors. 
We note that the increase of the ionic polarizability upon 
hydration in a high dielectric liquid has been previously 
observed in ab-initio calculations with PCM solvent [H[S]. 
This peculiarity was also revealed in Ref. [7] for water 
molecules, whose transfer from gas to liquid state was 
shown to be accompanied with a large amplification of 
their average dipole moment. In Section [HI B[ it will be 
shown that a similar hydration mechanism is present as 
well in ionic liquids. 

Moreover, in Fig. [3tc), one sees that the effective in- 
trinsic polarizability b V i exhibits in turn a monotonic de- 
crease upon hydration, until it reaches in the fully hy- 
drated state almost half of its gas phase value b P i . This 
indicates that upon hydration, the electronic cloud of the 



polarizable molecule increases in size, but also reaches an 
enhanced rigidity. In other words, the hydration opposes 
the electronic cloud deformation resulting from thermal 
fluctuations. Interestingly, comparison of Figs, k^ (a) and 
(c) shows that the increase of the electron cloud rigid- 
ity manifests itself at considerably lower concentrations 
than its expansion. Furthermore, in Figs. ^ (a) and (b), 
one notices that a significant departure from the MF be- 
havior with b m = and btot,i = b p i is observed above 
the characteristic solvent concentration p^ ~ 10~ 3 M. 
This shows that in Fig. pJFc) , the hardening of the elec- 
tronic cloud takes place already in the weak electrostatic 
coupling regime. Finally, in Figs. [3] (a)- (c), we note that 
although ions with a higher valency are clearly better 
solvated, the ionic dipole moment and polarizabilities ex- 
hibit weaker sensitivity to the molecular charge than the 
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FIG. 4: (Color online) Rescaled ionic dipole moment (dashed 
red curves) and total polarizability (solid black curves) 
against the gas phase polarizability b P i for the solvent den- 
sities (a) p s b = 2.0 x 10 -4 M and (b) p s b = 55.0 M, molecular 
charge |e»Ci| = 2, and solvent molecular size a — 1 A. The 
curves are from the full numerical calculation, the black and 
red circles respectively correspond to the limiting laws (231 
and ( |24[), and the black and red squares are from the expres- 
sions (|28| and ( 30 1 for the fully hydrated state. 



hydration energy in Eq. ( 13 ) characterized by a linear 



dependence on the charge [e^Cj | . 

In order to characterize the scaling of the hydrated 
polarizabilities with the gas phase polarizability b P i, we 
first consider the electrostatic weak coupling regime of 
dilute solvents. By evaluating in the dilute solvent 
regime the averages in Eqs. (20) and (22) at the order 
O y{n s a) 2 £B/a) , one obtains for the ionic dipole moment 
and the total polarizability 
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where we introduced the auxiliary functions 
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In Fig. [5| a), we compare the prediction of these 
asymptotic laws (circles) with the numerical solution of 
Eqs. (fl0|)-(fT3"|) (continuous curves) for a dilute liquid with 
density p s b = 2.0 x 10 -4 M. One notices that the behav- 
ior of the polarizabilities is characterized by two regimes 
separated by a peak located at b p i ~ a/3. Indeed, the 



asymptotic limit of Eqs. fl23| and (24) indicate that the 



average electronic cloud radius and total polarizability 
grow with the gas phase polarizability as b 
btot,i — b r 



bS and 



b i for b p i <C a/3 (left branch of the curves 



for 



in Fig. Uja)), and b mi ~ bj and b totli 
b pi ^> a/3 (right branch of the curves). Thus, the tran- 
sition between these two regimes results from a competi- 
tion between the solvent molecular size and the gas phase 
polarizability. 

In the opposite regime of concentrated solvents, the 
expansion of Eqs. ( 10 ) and ( 13 1 for K s a 3> 1 and b p /a <C 1 



yields for the hydration energy the asymptotic limit 
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(27) 



Neglecting the exponential term and expanding the to- 
tal distortion potential Ui(b) = hi(b) + ipi p (b) around the 
equilibrium position, we are left with the gaussian dis- 
tribution Ui(b) = (b—b m i) /(46^), with the average 
electronic cloud radius and effective intrinsic polarizabil- 
ity 
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Substituting these relations into Eq. (21), the total ionic 
polarizability follows as 
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Figur es [3{ a)-(c) show that the closed form expressions in 
Eqs. (28)-(|30l accurately reproduce the saturation val- 



ues of the ionic dipole moment and the polarizabilities 
(dashed horizontal curves). First of all, in Eq. (29), 



one notes that regardless of the ion charge, transferring 
the ion from the gaseous phase into the liquid environ- 
ment reduces its intrinsic polarizability by a factor three. 



Moreover, Eqs. (28) and (30) indicate that in the fully 
hydrated state, the ionic dipole moments and total polar- 
izability grow as the cubic root of the ion charge, which 
explains the weak dependence of the solvation on the 
molecular charge strength in Figs. |3][a)-(c). 

We compare in Fig. [l|b) the limiting laws ( 28 ) and ( 30 ) 
with the full numerical solution of the self-consistent 
equations for the solvent concentration p s i, = 55.0 M. 
These equations indicate that in the range b P i = 0.1 A to 
1.0 A, the dipole moment and polarizability of the fully 
hydrated ion grows with the gas phase polarizability ac- 

2 /3 

cording to the b ^ power law. We also note that interest- 



ingly, the hydrated polarizabilities in Eqs. (28)-(30) are 



independent of the solvent molecular size. This peculiar- 
ity stems from the fact that the complete hydration takes 
place in the parameter regime k s 3> a -1 , where the part 
of the dielectric susceptibility function associated with 
the rotation of solvent molecules (i.e. the third term on 
the r.h.s. of Eq. (Ill) makes no contribution to the hy- 
dration energy ipib(p} m Eq. (13). 



In our previous work on the MF theory of polar liq- 
uids at charged interfaces, it was shown that the non- 
local character of electrostatic interactions in the solvent 
results from the finite size of solvent molecules [20] • The 
effect of non-locality on the hydration mechanism can be 
estimated by varying the solvent molecular size a at fixed 
dipole moment p a = Qa. To this aim, we reexpress the 



dielectric permittivity function (11) in the form 



~<q) = 1 + 



JKsPo) 
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sm(qa) 
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(31) 



and calculate the total polarizabilities ( 20 )-( 22 ) with the 



above permittivity function by varying a with the dipole 
moment fixed at po = 1 A. In Figs, pi (a) and (b), the 
comparison of the curves with a = 1 A and 3 A shows 
that the increase of the solvent molecular size at fixed 
dipole moment lowers the average electronic cloud radius 
and the total ionic polarizability. Hence, non-locality 
weakens the hydration of the polarizable ion. To explain 
this peculiarity, we note that in the dilute ion regime, the 
inverse Fourier transform of the potential in Eq. ( 10 1 is 



given by a generalized Coulomb law, Vq(v) = £b/ [re(r 
with the local dielectric permittivity function 



e(r) = J 

z I JO 



dfc sin(fcr/a) 
T i(k) 



(32) 



and the adimensional wavevector k — qa. The dielectric 



permittivity profile of Eq. ( 32 ) is reported in Fig. k5 



First of all, it is seen that the close vicinity of the ion 
at r — is characterized by a dielectric void. Then, one 
notes that the dielectric permittivity function in Eq. ( 32 ) 



depends solely on the rescaled distance r/a. This means 
that an increase of the solvent molecular size amplifies 
the dielectric void around a polarizable molecule, and 



consequently reduces its hydration energy in Eq. ( 13 ) 



In the opposite point-dipole limit of solvent molecules 
a — > 0, the permittivity function (31) tends to the bulk 
permittivity, e (q) — > £& = 1 + AirlBPaPsb/S, and the hy- 
dration PMF (Il3| takes the simple form [2S] 



tpi P (b) = ip ip (b -> oo) + 



with the adimensional parameter 
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(34) 



Evaluating the integrals in Eq. (12 1 with the PMF (33 1, 



the moments of the electronic cloud oscillations can be 
expressed in terms of Meijer G-functions [27] . 
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The ionic dipole moment and total polarizability ob- 
tained from Eq. (35) is reported in Figs. [3la) and (b). 
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FIG. 5: (Color online) Effective dielectric permittivity of ionic 
liquids with bare polarizability b P i — 0.2 A and valency n. 
Solid curves are from the numerical solution of Eqs. ( |10|) -(|13|, 
dotted curves denote the full solvation limit in Eq. ( 46 1 , and 



(44 1 



square symbols from the approximative scheme of Eqs. 
(45 1. The black dashed curve is the MF dielectric permittivity 
Eb = 1 + £ P for n = 3. Inset : Total ionic polarizabilities (solid 
curves) and their saturation values from Eq. ( |30| ) (dashed hor- 
izontal curves). 



One notices that the point dipole result is very close 
to the case with finite solvent molecular size a = 1 A. 
Thus, for the model parameters chosen in this work, non- 
locality plays a minor role in the hydration process. It 
is interesting to note that in this parameter regime, the 
hydration of the polarizable ion can be solely described 
by the single coupling parameter T. 

By Taylor-expanding Eq. (35) in the regime r> 1, 
one obtains for the ionic dipol 



le moment and total polar- 
izability the following expressions, 



°pi = 2 r 1 / 3 + -r- 1 / 3 + o(r- 1 ) 

Opi o 
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6^6 
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(37) 



In Figs. f3^a) and (b), it is shown that the asymptotic 
laws p6]) and ( 37 ) can accurately reproduce the increase 



of the ionic dipole moment and total polarizability from 
p s b = 10 -3 M to complete hydration. These equations in- 
dicate that the fully hydrated state of the polarizable ion 
is reached with increasing solvent concentration through 
the gradual saturation of the parameter T in Eq. ( 34 ) . 



We consider next the counterpart of this hydration pro- 
cess in ionic liquids without solvent molecules. 



B. Cooperative solvation in ionic liquids 

Ionic liquids are promising salt solvents that gradually 
replace water in new generation energy storage devices 
such as graphene based capacitors |28) . The accurate 



knowledge of the dielectric permittivity of ionic liquids 
is needed to predict the charge storage ability of these 
devices. In ab-initio calculations of ionic liquids com- 
posed of small ions with negligible dipole moments [5], 
it was found that the contribution from electronic and 
orientational polarization of individual ions cannot alone 
explain the large dielectric permittivities measured in ex- 
periments 29 . Based on this observation, it was also 
argued that an additional polarization effect induced by 
the surrounding ions must be present to explain the high 
dielectric permittivity values. 

In order to shed light on this point, we consider in this 
part the closure equations (10)-(13) for an ionic liquid 



free of solvent molecules, and composed of two species 
of polarizable ions with the same bare polarizability b P i 
and bulk density pa,. Furthermore, the point charges on 
the polarizable molecules are e±.2 = ±1 and c\2 = ±n, 
which corresponds to the net molecular charges g x 2 = 
±(n— 1) (see Fig. [2th)). The dielectric permittivity of the 
medium at large separation distances from a central ion 
is obtained from the IR limit of Eq. (Ill, £& = e(q —¥ 0), 
and it is given by 



£b 
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where the total ionic polarizability defined in Eq. ( 22 ) has 
to be computed from the numerical solution of Eqs. ( 10 )- 



(13 1. Indeed, for an ionic liquid where the hydration of 



the polarizable ion affects the polarization of the sur- 
rounding medium in a self-consistent way, the solution of 
these equations is more tricky. Our numerical scheme 
consisted in solving these equations by iteration on a 
discretized Fourier lattice. Namely, at the first itera- 



tive level, the MF permittivity of Eq. ( 16 ) was used as 
the input function in the potential Eq. ( 10 ) in order to 



evaluate the hydration PMF in Eq. (13), and the latter 
was injected at the next step into Eq. (12) to obtain the 
updated dielectric permittivity function from Eq. (11). 



This procedure was continued until self-consistency was 
achieved. 

We illustrate in Fig. pj] the ionic polarizability (inset) 
and the dielectric permittivity of the liquid (main plot) 
obtained from the numerical solution of Eqs. (|To|-(13) 
(solid curves). First of all, it is seen that the increase of 
the ion density is accompanied with a strong amplifica- 
tion of the total ion polarizability, which in turn results in 
a rise of the dielectric permittivity of the medium. Then, 
in the inset of Fig. [5j we note that unlike the case of a po- 
larizable ion in a polar solvent (see Fig. [3ja)), the ionic 
polarizability and the full hydration density exhibits a 
pronounced dependence on the molecular charge. 

These effects can be shown to be driven by the self- 
consistent solvation of polarizable ions by their own field. 
To this aim, we introduce respectively the charge and 



dipolar screening parameters 



k 2 c = 4irl B J2 P^b1^ = 8Tr£ B (n - \f Plh (39) 

i 

K p = ^2 K % = ^^snptb, (40) 



and the corresponding coupling parameters £ c = {K c b P i) 
and £ p = (n p b p i) . In the dilute liquid regime, by expand- 
ing the closure relations ( 10 1 and ( 13 ) up to the order 



O (£ c ) and O (£ p ), one obtains for the solvation PMF 
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where we introduced the auxiliary function 
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One sees in Eq. (41) that the solvation energy is com- 



posed of a contribution from the charge screening (the 
first term on the r.h.s.), and a part resulting from the po- 
larizability induced dielectric screening of the ion by the 
surrounding ionic liquid (the second term on the r.h.s.). 
We display in Fig.[6]thc PMF of Eq. (41 ) for a monovalent 
ionic solution (n = 2) with concentration pn, = 5 x 10~ 5 . 
It is seen that in this dilute liquid regime, the charge 
and dielectric screening effects independently lower the 
bare distortion energy hi(b) with an equal weight, thus 
favoring the expansion of the electronic cloud. Then, we 
note that as in the case of a polarizable ion in a polar 
solvent considered in Section |IH A| the total distortion 
potential exhibits a minimum. In other words, in the 
liquid environment, the polarizable ion acquires a finite 
dipole moment. We emphasize that this effect has been 
previously observed in ab-initio calculations of ionic liq- 
uids composed of charges with fluctuating geometry [5] . 

In order to determine the relative weight of the di- 
electric and charge screening mechanisms in the renor- 
malization of the background dielectric permittivity be- 
yond the dilute regime, we will introduce an approxi- 
mative solution scheme of Eqs. ( 10 1-( 13 ). To this aim, 



we first redefine the hydration PMF of Eq. ( 13 1 by sub- 



tracting the constant energy in the dissociated state, 
Pipib) = ipip(b) — ipi P (b — > oo). Introducing the dimen- 
sioneless wavevector k = b P iq and separation distance 
x = b/b p i, this PMF can be expressed as 
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(43) 



where the statistical average of the functions inside the 
brackets is still evaluated according to Eq. (12) with the 
adimensional electronic cloud radius x' = b' /b p i as the 



10 



evaluating the PMF of Eq. (13) in this limit, we found 







FIG. 6: (Color online) Drude oscillator potential Eq. ([T]) (blue 
curve), the first (charge screening) and second term (dielec- 
tric screening) of the ionic solvation energy in Eq. ( |41[ ) de- 
noted respectively by the red dotted and dashed curves, and 
the total distortion potential (black curve) for an ionic liq- 
uid composed of polarizable ions only, with ionic density per 
species ptb = 5 x 10 -5 M, gas phase polarizability b V i — 1 A, 
and molecular charge |ejCj| =2. 



integration variable. We now assume that the hydration 
PMF affects the electron cloud mainly at small separa- 
tions x < 1. This implies that in Eq. (43), only small 
wavevectors k < 1 make a significant contribution to 
the integral. Based on this assumption, by expanding 
the sinusodidal functions inside the bracket of Eq. ( 43 ) 
at the order 0(fc 2 ), the integral can be evaluated ex- 
actly. Within this approximation, the complicated inte- 
gral equations (fl0|-([l3| for the dielectric permittivity are 



reduced to a simpler non-linear equation, 



£6 = 1 



Z p f™dxx 4 e- x2 / 4 -^^ 
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where we made use of Eqs. (12) and (38). 

In Fig. [5j it is shown that the numerical solution of 



Eq. (44) can accurately reproduce the dielectric permit- 



Dtained from the closure equations (10)-((T3|) over 



tivity o 

the whole density range. We now note that in the solva- 
tion PMF of Eq. ( 45 ) , the contribution from the dielectric 



and charge screenings correspond respectively to the first 
constant term £{, and the second exponential function in 
the numerator. This equation indicates that while in- 
creasing the ion concentration from the dilute regime, 
the exponential term is gradually dominated by the con- 
stant term in the numerator and becomes negligible for 
£b 3> 1. Thus, charge screening makes a significant con- 
tribution to the dielectric permittivity exclusively at low 
ion concentrations. 

To asccertain the latter point, we now consider the 
strict limit of large liquid densities with K p b p i 3> 1. By 



that the total ionic polarizability is still given by the 
expression (30) (see the horizontal lines in the inset of 
Fig. [5]). Substituting this relation into Eq. p8| , one ob- 
tains the dielectric permittivity of the ionic liquid at the 
fully solvated state 
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In the main plot of Fig.[5j it is shown that this closed form 
expression is a very good approximation for the dielectric 
permittivity of the ionic liquid beyond the dilute regime. 
One can note that in Eq. (46), the dependence of the 
permittivity on the charge screening parameter £ c has 
disappeared. This shows that close to the full solvation 
state, the collective solvation mechanism is solely driven 
by the dielectric screening induced by polarizable ions. 

We also compare in Fig. [5] the MF level bulk dielectric 
permittivity e\, = l+£ p for the ion valency n = 3 with the 
self-consistent result. The MF theory that neglects the 
collective ionic solvation is shown to strongly underesti- 
mate the dielectric permittivity of the ionic liquid. This 
observation is in line with Ref. [2] where the rotational 
polarizability associated with the gas phase dipole mo- 
ment of ions was shown to be unsufficient to explain the 
high dielectric permittivity of ionic liquids. This suggests 
that the cooperative hydration mechanism scrutinized in 
this part brings the main contribution to the dielectric 
permittivities of ionic liquids. Hence, correlation effects 
cannot be neglected in polarizable liquids. 



IV. CONCLUSION 

We have introduced in this article a classical electro- 
static theory of polarizable ions in high dielectric liq- 
uids. Within this theoretical framework, we have scru- 
tinized the physical mechanism behind the ionic solva- 
tion properties observed in ab-initio calculations of polar 
solvents [H [5] and ionic liquids [HE]. In the first part 
of the article, we presented the electrostatic formulation 
of polarizable ions immersed in polar solvents composed 
of dipolar molecules with finite size. Then, we derived 
from the Dyson equation the electrostatic self-consistent 
relations accounting for the electrostatic correlations be- 
tween the particles in the liquid. 

The second part of the article was devoted to the hy- 
dration of a single polarizable in a polar solvent such as 
water. It was shown that the electrostatic energy release 
experienced by the polarizable ion upon hydration re- 
sults in the expansion of its electronic cloud. As a result, 
the ion carrying zero dipole moment in the gas phase 
acquires in the liquid environment an average dipole mo- 
ment. However, the hydration also amplifies the rigidity 
of the electronic cloud, thereby opposing its deformation 
induced by thermal fluctuations. In qualitative agree- 
ment with quantum molecular calculations with PCM 
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solvent jH [5] , the overall effect was shown to be an en- 
hancement of the gas phase polarizability upon hydra- 
tion. 

In the third part of the article, we have investigated 
a cooperative solvation mechanism in ionic liquids free 
of solvent molecules. We have found that similar to the 
case of a polarizable ion in the polar solvent and in agree- 
ment with ab-initio calculations of ionic liquids [8 , each 
polarizable ion acquires in the liquid a finite dipole mo- 
ment and an increased polarizability. This effect resulting 
from the polarization field generated by the surrounding 
ions self-consistently amplifies the dielectric permittiv- 
ity of the medium. We note that this solvation induced 
amplification of the dielectric permittivity is substantial 
even in the weak electrostatic coupling regime of dilute 
liquids. This suggests that the self-consistent solvation 
mechanism brings the dominant contribution to the di- 
electric permittivity of ionic liquids composed of small 
ions with negligible permanent dipole moment in the gas 
phase P]. 

We have introduced the first microscopic theory of 
ionic hydration in explicit solvent, and we emphasize 
that the model as well as the theoretical scheme need 
refinements. First of all, it should be noted that our 
approach does not account for the hydrogen bond for- 
mation in water solvent, which is believed to amplify the 
dielectric permittivity of water |30| . This complication 
expected to become significant beyond the dilute liquid 
regime should be addressed in a future work by extend- 
ing our approach beyond the gaussian field approxima- 
tion, i.e. by opting for a more sophisticated closure to 
solve Eq. ^. An additional complication for solvents at 
physiological concentrations comes from the importance 
of excluded volume effects associated with the finite size 
of the particles in the liquid. The first step to generalize 
the model in this direction consists in including simple 
hard-core or repulsive Yukawa interactions between the 
particles as in Refs. [3TH55] . Then, our theoretical scheme 
should be extended to a second order cumulant expan- 
sion of the grand potential around the reference Hamil- 



tonian Eq. ^ . This generalization would allow to deter- 
mine how much our results are quantitatively modified 
beyond the dilute liquid regime. Indeed, we expect hard- 
core interactions between solvent molecules and ions to 
reduce the polarizability increase induced by the electro- 
static hydration mechanism. In this sense, the results 
presented in this article beyond the dilute solvent regime 
should be considered as an upper boundary for the ac- 
tual ionic cloud expansion effect. Our results should be 
also compared at the next step with MC simulations of 
the polarizable ion model introduced in Sec.[TTJ but these 
simulations are currently unavailable. 

Finally, the consideration of the induced polarizability 
with a classical Drude potential is another limitation of 
the present model. Actually, it should be noted that 
the ionic dipole moments in the solvated state provided 
by our theory are larger than the values observed in 
ab-initio calculations [U |8]. For example, the Pauli 
exclusion effect neglected by the classical approach is 
expected to partially suppress the hydration induced 
expansion of the electron cloud. However, refinements at 
the quantum level are of course beyond the scope and the 
main message of the present work. Indeed, the ability 
of the theory to qualitatively capture ionic hydration 
effects observed in quantum molecular calculations for 
both polar solvents and ionic liquids on the one hand, 
and the presence of these effects in the dilute liquid 
regime where the complications discussed so far are 
not expected on the other hand confirm the physical 
consistency of the model with real dielectric liquids. 
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